The eccentricities of flow - elliptic flow fluctuations and evidence for transverse 
localization in the initial state of the matter in relativistic heavy ion collisions* 
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Recent measurements of event-by-event elliptic flow in Au+Au coUisions at -^snn ~ 200 GeV 
exhibit large relative fluctuations of about 40-50%. The data are well described by fluctuations in 
the shape of the initial collision region, as estimated event-by-event with the participant eccentricity 
using Glauber Monte Carlo. These results, combined with the demonstrated participant eccentricity 
scaling of the elliptic flow across nuclear species, constitute evidence of transverse granularity in the 
initial matter production in these collisions. 



I. INTRODUCTION 



Elliptic flow {V2) is one of the key observables in the un- 
derstanding of the dynamics of heavy ion collisions. The 
observation of a significant azimuthal anisotropy in the 
momentum and/or spatial distributions of the detected 
particles relative to the reaction plane, is direct evidence 
of interactions between the initially produced particles 
in heavy ion collisions. These interactions must occur 
at relatively early times, since expansion of the source 



* Presented at the International Workshop on Hadron Physics and 
Properties of High Baryon Density Matter, Xi'an, China, November 
2006. 



rapidly reduces the magnitude of the spatial asymmetry. 

Typically, the connection between the initial and final- 
state anisotropy is provided by hydrodynamical models 
that relate a given, initial source shape to the distribution 
of produced particles. In such calculations, it is common 
to use smooth, event-averaged, initial conditions. How- 
ever, event-by-event fluctuations in the shape of initial 
interaction region must not be neglected. As a means 
to quantify the effect of initial-state eccentricity fluctu- 
ations, PHOBOS has introduced the "participant eccen- 
tricity" The magnitude and shape of epart as a 
function of centrality were found to be robust to varia- 
tions of the Glauber parameters. 

Fluctuations in the shape of the initial state interac- 
tion region might be expected to be more pronounced 
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for smaller species and for more peripheral events. This 
is borne out by the data, where a striking agreement 
between the elliptic flow signals in Cu+Cu and Au+Au 
data as a function of centrality is obtained when scaled 
by the participant eccentricity, V2/{€part) [llj Si- 
Given the obvious importance of fluctuations in the ini- 
tial state to the final elliptic flow measurement, PHOBOS 
recently developed a new technique to measure dynami- 
cal flow fluctuations in our data 0, . The new analysis 
and results stemming from it are summarized below. 



II. ELLIPTIC FLOW FLUCTUATIONS 
ANALYSIS TECHNIQUE AND RESULTS 

The PHOBOS flow fluctuations analysis extracts the 
flow signal, event-by-event, from data in the PHOBOS 
multiplicity array which detects a very large fraction 
of the produced particles over the pseudorapidity range 
1 77 1 < 5.4 Q. The event-by-event measurement is done 
using a maximum likelihood fit with two parameters to 
the hit information over the full acceptance of pseudo- 
rapidity. The parameters determined in the maximum 
likelihood fit are the observed elliptic fiow at midra- 



distribution of V2, g(w2 ); by 



pidity, W2(0), and the reaction plane angle. 
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pseudorapidity dependence in the likelihood fit, V2{ri), is 
parametrized with a triangular shape, V2"{'r]) = V2 (1—-'^), 
or alternatively with a trapezoidal shape, v^^'^^lrj) = 
{f4''''(l)t?f|r,|>2 ' "^here V2 = W2(0). Both parametriza- 
tions yield good descriptions of the previously mea- 
sured (mean) V2 {r]) shapes [Q] . The fit includes a proba- 
bility density function that corrects for non-uniformities 
in the acceptance of the used sub-detectors. 

In order to disentangle known (mostly statistical) from 
unknown (dynamical) contributions to the measured fiow 
fluctuations, a detailed knowledge of the detector re- 
sponse is required. A response function, K{v2^^ ,V2,n), 
is defined as the distribution of the event-by-event ob- 
served elliptic flow, V2^^, for events with constant input 
flow value, V2, and multiplicity, n. This response function 
is designed to account for detector deflciencies, as well 
as for multiplicity and finite-number fluctuations when 
mapping the true distribution of V2 to the distribution 
of ^2'^'^. If f{v2) is the true V2 distribution for a set of 
events in a given centrality class, f{v2) is related to the 
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K{vf^V2,n)f{v2)N{n)dv2dn, (1) 



where N(n) is the multiplicity distribution of the events 
in the given set of events. 

To obtain the kernel in bins of 112 and n with enough 
precision, would require on the order of 100 million MC 
events. Instead the kernel can be parametrized, allowing 
the use of about 1.5% of that statistics to reach the re- 
quired precision. For a perfect detector, the response is 
given by Eq. (A13) from Ref. Q (with a^vf''', a^V2 
and AI n). In practice, however, it turns out that V2 
is suppressed, with the suppression dependent on n, and 
that the resolution (a) has a constant background con- 
tribution. With ^2"^ = {An + B)v2 and a = C/^/n + D, 
this leads to 



K{v^'^^,V2,n) = ^exp[- 
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(2) 



where Iq is a modified Bessel function. The four unknown 
parameters {A,B,C,D) are obtained using the modified 
HIJING samples. 

In order to determine the mean and variance of the 
true V2 distribution, f{v2), and extract the fluctuations, 
we assume a Gaussian distribution for f{v2), with two 
parameters, (1)2) and ay.^. For given values of {V2) and 
7 it is possible to take the integral in Equation ([1]) 
to obtain the expected distribution, gcxp{v2^^\{v2) ,(Jv2) ■ 
By comparing the expected and observed distributions, 
the values for {V2) and found by a maximum 

likelihood flt. 

The analysis chain is applied to the ^snn — 200 GeV 
Au+Au data set from Run 4 in bins of centrality. The 
results {{V2) and cr^^) are obtained separately for trian- 
gular and trapezoidal V2{r]) shape and averaged over 10 
bins of collision vertex (2 cm width). The systematic 
errors (all sources added in quadrature) are estimated by 
including variations from different vertex and bins 
and changes introduced by the triangular and trapezoidal 
V2{ri) shapes. Since the functional form of the true dis- 
tribution is unknown, also differences arising from a flat 
rather than a Gaussian ansatz for /(W2) are included. 
Furthermore, we have performed extensive studies of the 
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analysis response to MC samples, prepared to match 
the observed multiplicity distribution and (^2), in bins 
of centrality and known input values of fluctuations (in- 
cluding zero) . The discrepancies of the obtained to the 
input fluctuations in the vicinity (±20%) of the mea- 
sured values per centrality bin are an additional source 
of the systematic error. Their contribution becomes in- 
creasingly important for lower values of V2, which is the 
dominant reason to exclude the 0-6% most central bin 
from the analysis. Fig.[T]shows the mean flow results ob- 
tained with the event-by-event analysis projected to mid- 
rapidity by {V2) — 0.5 (yi (U2") + (^2'''''))' which agree 
well, within scale errors, with the published results of 
the event-averaged sub-event based technique using hits 
or tracks This constitutes an important, completely 
independent verification of the event-by-event analysis. 
Fig. [2] presents the relative flow fluctuations, 0-^2/ {v2)- 
In the ratio most of the forementioned systematic errors 
scale out. Large relative fluctuations of about 40-50% 
are observed, with relative little centrality dependence. 
The PHOBOS results presented here are consistent with 
a recent, preliminary measurement of elliptic flow fluctu- 
ations by the STAR collaboration 0] . 

III. FLUCTUATIONS EXPECTED FROM THE 
PARTICIPANT ECCENTRICITY MODEL 



pant eccentricity obtained in MCG simulations (to deter- 
mine the 90% confidence level band shown in Fig. [21 the 
Glauber parameters were varied within reasonable limits 
as described in Ref. Q). The contribution from iVpart 
fiuctuations, estimated using a fit to the (U2) data and 
the known iVpait distributions from the PHOBOS central- 
ity trigger studies, is neglible in the measured centrality 
range. 



CM 



0.1 

0.08 
0.06 
0.04 
0.02 



PHOBOS preliminary 

Au+Au, 200 GeV, If, |ri|<1 



\ ^ — 



• Event-by-event 
a Hit-based 
H Track-based 



50 100 150 200 250 300 350 



N 



part 



FIG. 1: Mean flow, (U2), as a function of centrality, for 
y/swN = 200 GeV collisions at mid-rapidity, measured by the 
event-by-event analysis, compared to the published results 
obtained with the event-averaged, sub-event based technique 
using hits or tracks 



The participant eccentricity picture accounts for 
nucleon-position fiuctuations in the participating nucleon 
distributions by calculating the eccentricity, event-by- 
event, with respect to the principal axes of the overlap 
ellipse in a MC Glauber (MCG) simulation. In a hy- 
drodynamical scenario, such fluctuations in the shape of 
the initial collision region would lead naturally to cor- 
responding fluctuations in the elliptic flow signal. To 
estimate their magnitude, it is assumed that V2 0ce event- 
by-event. This leads to a^^/ {^2) = fe/ (e), where a^^ (ce) 
is the standard deviation of the event-by-event distribu- 
tion V2 (e) , provided there are no other sources of elliptic 
flow fluctuations. Neglecting all other sources of elliptic 
flow fluctuations, the participant eccentricity MCG simu- 
lation predicts relative fluctuations {<^v2l{'^2)) of 35-50% 
in Au+Au collisions at ^snn = 200 GeV. The prediction 
is shown in Fig. [2] as a function of A'pait- The data are 
well described by the prediction, a^/ (e), from the partici- 



0.8 
0.6 
c^ 0.4 
0.2 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
PHOBOS preliminary 




Au+Au, 200 GeV, If, |ti|<1 


- \ 

\ 




\ 


• 


■I » • • • 


• " X — 


~_ • Data (90% G.L.) 




- ^ MCG o/(G)(part) 


^^^_0„/(v>(N ) " 

^2 2 part 
1 , , , , 1 , , , , 1 , , , , 1 , , , 



% 50 100 150 200 250 300 350 
'^^'part 

FIG. 2: Relative flow fluctuations, 0-^2/ {""2), as a function 
of centrality, for ^snn = 200 GeV collisions at mid-rapidity, 
compared to the prediction, a^/ (e), from the participant ec- 
centricity, and to an estimate of A^'part induced fluctuations 
using a fit of {V2} (A^part) 
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IV. CONCLUSIONS 

Recent results on event-by-event elliptic flow fluctua- 
tions in Au+ Au collisions at ^/sfm = 200 GeV exhibit 
large relative fluctuations of about 40-50%, relatively in- 
dependent of centrality. The new data are well described 
by fluctuations in the shape of the initial collision region, 
as predicted with the participant eccentricity using MCG 
simulations. 

These results substantiate conclusions from previous 
studies by PHOBOS on the relevance of such evcnt- 
by-event fluctuations for the elliptic flow across nuclear 
species. The initial-state geometry seems to drive the 
hydrodynamic evolution of the system, not only on aver- 
age, but event-by-event. The success of the participant 
eccentricity model in describing both the geometric scal- 



ing and fluctuations of the elliptic flow is consistent with 
the matter present in the initial stage of relativistic heavy 
ion collisions being created with a transverse granularity 
similar to that of the participating nucleons. 
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